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Fluctuation relations are derived in systems where the spin degree of freedom and magnetic interactions 
play a crucial role. The form of the non-equilibrium fluctuation theorems relies in the assumption of a local 
balance condition. We demonstrate that in some cases the presence of magnetic interactions violates this con- 
dition. Nevertheless, fluctuation relations can be obtained from the micro-reversibility principle sustained only 
at equilibrium as a symmetry of the cumulant generating function for spin currents. We illustrate the spintronic 
fluctuation relations for a quantum dot coupled to partially polarized helical edges states . 
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Introduction. Non-equilibrium fluctuation theorems 
(FTs) 01-0], widely used for macroscopic systems, are based 
on the thermodynamics governing the physical processes 
when they evolve forward and backward in time. The bound- 
ary conditions for the forward and the time-reversed proces ses 
determine the balance condition for the entropy exchange 
and therefore the form of the fluctuation theorem Q3J] . The 
applicability of the non-equilibrium FTs to quantum systems 
has become an exciting problem and, in particular, to the case 
of the charge transfer phenomena in mesoscopic systems in 
the context of the full counting statistics Interestingly, 
relations akin to the fluctuation-dissipation theorem 10-11] 
have been formulated beyond the linear response regime 
10, [HUH] • These fluctuation relations relate nonequilibrium 
fluctuation and dissipation coefficients for phase-coherent 
conductors. However, the role of a genuine quantum property 
such as the spin degree of freedom in the fluctuation relations 
has not been yet investigated in detail. Our motivation is not 
only fundamental since the electronic spin offers enormous 
advantages to create devices with unusual and extraordinary 
new functionalities [23]. The purpose of this work is thus to 
generalize the fluctuation relations for spintronic systems. 

Fluctuation relations are generated from the cumulant gen- 
erating function (CGF) F{\) = In £ Q P(Q, t)e~ ixQ where 
P(Q) is the charge distribution function. Firstly, the CGF 
T is expanded in a Taylor expansion in terms of affinities 
A = {qVi/ksT^ qV2/kBT : • • • ) (q is the electron charge, 
ks is the Boltzmann constant, T is the temperature and Vi, 
i = 1,2 .. . are the applied voltages) and counting fields 
X = (xij X2: • • • ) around the equilibrium condition. Then, 
thanks to the symmetries F(0,A) — (probability conserva- 
tion condition) and F(—A,A) = (global detailed balance 
condition) fluctuation relations among the higher-order non- 
linear cumulants are found. Indeed, the symmetries of T can 
be considered as the non-equilibrium FT versions for the cur- 
rents within a transport theory. Initial experiments by using a 
mesoscopic dot interferometer have tested these relations 112411 . 
In this experiment, the noise susceptibility and the second or- 
der conductance were found to be proportionally related. 

Spins are sensitive to magnetic fields and also to electric 
fields due to spin-orbit interactions. Fluctuation relations for 
the charge transport have been formulated in the presence of 
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FIG. 1. (Color online) (a) Sketch of a quasi-localized level (s a ) cou- 
pled to chiral edge states with Ti,2 and driven out of equilibrium with 
Vl and Vr bias voltages. For B > the upper (lower) edge state is 
injected from Vl (Vr). (b) Localized level coupled to unequally spin 
populated helical edge states due to the spin injection from the ferro- 
magnetic leads. Then, the tunnelling couplings are spin-dependent: 
r(i,2)t(i)- Ferromagnetic electrodes: larger light area (green) cor- 
responds to majority spins whereas smaller dark area (purple) are 
minority spins. Helical states: upper left (upper right) movers are 
spin up (spin down) carriers injected from Vr (Vl). 



magnetic fields Ill4 -|16ll. In Ref. 11411 the non-equilibrium 
FT for the forward and backward charge distribution proba- 
bility at opposite B polarities P(Q, B)/P(-Q, -B) = e QA 
was used to derive such fluctuation relations. However, some 
caution is needed since P(Q,B) and P(—Q,—B) are con- 
sidered for a system driven out of equilibrium in which the 
interacting internal potentials are no longer even functions of 
B |25 1 and the application of such theorem may break down 
fl5n . To circumvent this obstacle, Ref. lfl5n uses a symme- 
try of T associated with the micro-reversibility condition only 
at equilibrium P(Q, B)a^o = P(— Q, — P)a^o- We here 
derive the spintronic fluctuation relations in the same spirit 
when time-reversal symmetry is broken not only by exter- 



nal magnetic fields but also by the presence of ferromagnetic 
electrodes. In this case, at equilibrium P(Q, B, p)a^q = 
P(Q,—B, ~p)a^o, where p is the lead magnetization [26]. 

We illustrate our findings with a quasi-localized level cou- 
pled to helical edge states which are partially polarized by 
the presence of polarized electrodes [see Fig. 1 (b)]. Helical 
modes have been observed in topological insulators [27] and 
proposed to occur in quantum wires [28] and in carbon nan- 
otubes ll29ll . This quantum spin Hall state consists of gapless 
excitations that exist at the boundaries in which its propaga- 
tion direction is correlated with its spin due to the spin-orbit 
interaction. By electrostatic gating, quasi-localized states can 
form in the interior of the carbon nanotubes and quantum 
wires. Furthermore, ferromagnetic contacts have been suc- 
cessfully attached to these nanodevices [0. Finally, Ref . M\ 
suggests to create a quasi-bound state in quantum spin Hall se- 
tups by using ferromagnetic insulators that serve as tunnelling 
barriers. 

Local Detailed Balance. Consider a system described by 
a set of m discrete states coupled to ^-electronic reservoirs. 
We assume that its dynamics is governed by the master equa- 
tion dp/dt = Wp, where W is the transition rate matrix, 
and p denotes the occupation probabilities for the m- states. 
The exchange of energy (AEW) or particles (AN&) in the 
^-th-reservoir with inverse temperature is described by 
adding counting fields (x E , Xn) t0 me off-diagonal matrix 
elements of W. Thus, for the upper off-diagonal the tran- 
sition rate from the state n to the state m are modified ac- 
cording to W nm = W^e** aeW anW (for n < 
m) whereas for the lower off-diagonal terms these rates are 



W nrn = YiiW^e-XB^'-xX**" (n > m ). Usually, 
boundary conditions are taken into account through the lo- 
cal detailed balance (LDB) condition in which weight factors 
e -^ l) CHt-nNt) (<ft £9 and Nl denote the Hamiltonian and the 
particle number operator, respectively for the ^th-reservoir) 
balance forward and backward processes. To be more spe- 
cific: 
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From the LDB condition the equality W(xe\x$) = 
W T (f3^ - Xe\ -P w » w - Xn ] ) is automatically satisfied 
reflecting the following symmetry for the generating function 
T [which is constructed from W(xe > Xtv*)! 1 
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Although in many systems we can assume some type of 
LDB condition, in general, Eq.(Q} is not fulfilled [32]. To 
see this in a quantum conductor, we consider the system 
sketched in Fig. 1(a) in which the presence of a magnetic 
field breaks time-reversal symmetry. The system consists of 
a quasi-localized state with energy Sd in the Coulomb block- 
ade regime coupled to two chiral states propagating along the 
opposite edges of a quantum Hall conductor (filling factor 




FIG. 2. (Color online) Electrostatic model for the quasi-localized 
state connected to helical spin edge states. Spin up(down) localized 
level is capacitatively coupled to the upper up(down) helical channel 
with capacitances C u i( U 3)[Cdi(da)] an d to the lower up(down) he- 
lical channel with capacitances C u2 ( U 4)[Cd2(d4)]- A mutual capac- 
itance between up and down localized levels is accounted for with 
C. 0t(4O denotes the spin up (down) internal potential for the quasi- 
bound state. 



v = 1) |[T3, 0, [33|]. In the infinite charging energy limit 
case only two dot charge states are permitted: |0) and |1). For 
positive magnetic fields B > carriers in the upper (lower) 
edge state move from the left (right) terminal to the right (left) 
terminal. The current flow is reversed for B < 0. Inter- 
action between the quasi-localized state and the edge states 
takes place via tunnel couplings F\ and T2 and capacitive 
couplings C\ and C2. The chiral coupling involves differ- 
ent transition rates depending on the polarity of the magnetic 



field. For a positive B we have W ( 
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Ti(2)[l - f(B,fjL L ( R) )], where f(B,p L{R) ) 



1 / (1 +exp (3[pd(B) — Pl(r)]) is me Fermi-Dirac distribution 
function, /iz^i?) = qVum + Ep denotes the electrochemical 
potential in the lead L(R) with Ep = the Fermi energy, and 
Pd is the electrochemical potential of the quasi-localized state 
which is self-consistently calculated and depends on the B 
orientation |0]. ¥ov B > 0,p d (B)-p L = e d -(l-rj)qV/2, 
where r\ = {C\ — C<i)l{C\ + C2) is the capacitance asym- 
metry parameter and V = Vl — Vr. For B < the mo- 
tion of the edge states is reversed and then pd(—B) — pl = 
ed-(l + rj)qV/2. Because of the fact that p d ( B) ^ p d (~B), 
the LDB condition is not satisfied. Clearly, 
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where /? is the common inverse temperature and f eq = 1 — 
2 feq with f eq being the Fermi function at equilibrium (Vl = 
Vr) . Importantly, the violation of the LDB condition occurs 
for asymmetric capacitance couplings only. In the symmetric 
case or at equilibrium, Eq. (1) is recovered. The violations of 
LDB are thus a consequence of asymmetric, chiral states out 
of equilirium 

We now show that violations of the LDB conditions are 
also present in the absence of magnetic fields and when the 
spin degree of freedom is explicitly accounted for. For that 
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purpose we consider the system sketched in Fig. 1(b), a quasi- 
bound state which is tunnel coupled to helical edge states. The 
helical modes are partially polarized due to their coupling to 
two ferromagnetic electrodes with parallel magnetization and 
equal polarization p. In this manner, polarized helical edge 
states are described with a spin-dependent density of states 
(DOS) D ia = (1 + sp)Di/2, where s = +(-) for t (I)- 
helical mode and, Di, with i = u,d denoting the upper and 
lower edge-state DOS in the absence of polarization 0411 . In 
the wide band limit approximation, the tunnelling rates be- 
come spin-dependent, T is = 7r\ti\ 2 D is , with \ti\ 2 the tunnel 
probability from the i-th edge state. Defining I\ = 7r\ti\ 2 Di, 
wefindT is = (l + sp)Ti/2. 

Our transport description also includes an electrostatic 
model for interactions between the dot and the edge states. 
Within the mean-field approach, the electrochemical capaci- 
tive coupling C M consists of a geometrical capacitance contri- 
bution, C g , which depends on the width and the height of the 
tunnel barrier, and a quantum capacitance term, Ci S , q , which 
we take as proportional to D is : 
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We emphasize that the capacitive couplings are, in general, 
spin dependent l35ll . For sufficiently large geometrical capac- 
itances, C g ^> Ci S ^ q , we find from Eq. (4) the capacitative 
couplings 

1 +p 1 — p 
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where d = q 2 Di with C ul ( dl ), C u2 ( d 2) being the capaci- 
tive couplings between left (right) movers with up (down) spin 
along the top edge and the dot electron with spin | (|) whereas 
C U 3(d3) an d C U 4(d4) couple the same dot state with right (left) 
movers along the bottom edge with up (down) spins (see Fig. 
2). For thin edge states, Di depends on the steep confine- 
ment potential at the top and bottom edges, which will gen- 
erally differ Q. Then, we take d = C 3 and C 2 = C 4 
but C\ 7^ C2. As a consequence, the capacitive couplings 
between the dot and the edges is asymmetric: r\ ^ 0. Fur- 
thermore, since the upper and lower edge modes are equally 
polarized, one has Cdj — C u j . 

Consider for the moment the case where the capacitance 
coupling between the dot states is neglected (C = 0). Then, 
we calculate the spin- dependent electrochemical potential of 
the dot and find the simple relation fi a (p) — /^a(-p) = prjV. 
Now, in a time-reversal operation we have to invert the lead 
polarization, the edge state spin index, and the dot spin. Doing 
so, we obtain an invariant result only at equilibrium ( V = 0) 
or for symmetric capacitive couplings. But, in general, when 
V ^ the original state is not restored and, as a consequence, 
LDB is not fulfilled: 

= e*"-^[l - Pvpffe, + C(V% (6) 

where a = {t, 1} denotes the dot spin index. We stress that 
helicity is needed in our example to find departures from LDB . 
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Although we cannot rule out the possibility that nonchiral, 
spintronic systems (e.g., a dot directly attached to ferromag- 
netic leads) might show such departures if coherent tunneling 
or strong correlations are taken into accoutn, our conceptually 
simple system already exhibits the effect with fully analytical 
expressions. 

Spintronic fluctuation relations. We now treat on equal 
footing the presence of both magnetic fields and polarized 
contacts. The spin-dependent probability distribution satisfies 
the micro-reversibility condition, but only at equilibrium 
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where a and s are the lead and spin indices, respectively and 
p = (PjP'j---) contains the magnetizations for the leads. The 
CGF ^({ix}-) A) can be expanded in terms of powers of volt- 
ages and counting fields 
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f {kas}t{ i a }=\{d k - +la H{ix},A)/d{ix as ) k -dA l a 

as 

(9) 

where k and I are non-negative integers. From the deriva- 
tives of J r ({ix}j A) with respect to the counting fields, the 
cumulants are generated. In this way, the average current 
through terminal a with spin s is derived from Eq. ©as 
(las) — f{i as }- Similarly, second cumulant (current-current 
correlation) S as ^ s > = (AI a3 AIp 8 *) (AI as = I as - (I a8 )) 9 
where I denotes the current operator) and the third cummu- 
lant C as p s > ls » = (AI as AIp s f AI 7S ff) are given by S asps' = 
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expand, both, the current (I as ), and the noise (S a3 p 3 /) in 
powers of the applied voltages as follows 
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Here G 
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ductance, G { ^ s jk 



f{i as },{ij} corresponds to the linear con- 
f{i as },{iji k } is the second-order con- 
ductance, and S^lp^j = f{i as i ps ,},{i 3 } is the noise sus- 
ceptibility. Fluctuation relations are expressions that relate 
the /-coefficients at different order in voltage. To derive 
explicitly these relations we employ the micro-reversibility 
condition at equilibrium F(iXa S , i>Xps'i ' ' ' > A +^)U=o = 
F(-iXas, -iXpZ>, "' ,A -B)\ A =o [cf. Eq. ©]. It is con- 
venient to define the symmetrized (+) and anti- symmetrized 
(— ) combination of the /-factors 



ffk. 



{kas} ,{lj} 



f{k a *},{h}( B >P) ± f{k a s},{lj}(~ B > 



-v). 

(11) 



4 



where f{k aS },{ij}(~B, —p) is generated by means of time re- 
versal operation B — > —B,p^ —p, and s — > —8. According 
to Eq. d7]) the / ± -factors are even(odd) functions under time- 
reversal operation. This even-odd property is translated into 
the following relations for the equilibrium coefficients [in the 
sense of a voltage expansion, see Eqs. dTOtl: 
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Now by using the global detailed balance condition 
J r (—A,A)± = 0, and the probability conservation 
J 7 (0,A)± = one can derive the spintronic fluctuation re- 
lations among different f± -factors. Here we explicitly show 
those that relate the coefficients appearing in the third cumu- 
lant, noise and the conductances in the voltage expansion of 
Eq. GU: 
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Fluctuation relations between even higher-order response co- 
efficients toward the strongly nonequilibrium domain can be 
similarly found, relating different current cumulants at differ- 
ent order; however, the resulting expressions, already in the 
spinless case, look rather cumbersome 111511 . 

We verify Eq. O in a multi-terminal setup in which 
LDB condition is broken. For that purpose we general- 
ize the two terminal quantum spin Hall bar system [Fig. 
1(b)] to the multi-terminal case in which upper and lower 
helical modes are now connected to different terminals Vi, 
i = 1 • • - 4 [see inset in Fig. 3(d)]. We additionally con- 
sider spin-flip relaxation events within the quasi-bound state 
that can occur due to spin- spin interactions with a spin fluc- 
tuating environment (hyperfine interaction, spin-orbit inter- 
actions, etc.). We phenomenologically model this rate as 
7jJ = 7s/ exp [(e a — £&) / '(2/^bT)]. Notice that due to spin- 
flip events spin up and down currents are correlated and then 
Eq. (fT3l) is satisfied in a non-trivial manner. We emphasize 
that Eq. (fT3l) is verified (see Fig. 3) even for a finite capaci- 
tance asymmetry where the LDB condition is not met. 

Conclusions. In short, we have shown that the applica- 
bility of non-equilibrium FT when magnetic interactions are 
present is not a priori ensured. We illustrate this state- 
ment by using a quasi-localized level coupled to a chiral one- 
dimensional conducting channels. We demonstrate that lo- 
cal detailed balance condition is not satisfied when a mag- 
netic field is included and the system is driven out of equilib- 
rium. Importantly, we have derived the fluctuation relations 
for spintronic systems and have explicitly verified them in the 
illustrative case of a quasi-localized state coupled to partially 
polarized helical edge states. Our formalism is based on zero- 
frequency fluctuations and time-independent fields but in the 
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FIG. 3. Verification of spintronic fluctuation relations as a function of 
7 S / in the presence of magnetic interactions, B, and p for different 
values of polarization p: (a) p = 0, (b) p = 0.25, (c) p = 0.5, 
and (d) p = 0.75. Upper helical modes: left (right) movers are 
spin up (down) injected with voltage Vi( 3 ). Lower helical modes: 
right (left) movers are spin up (down) injected with voltage V 2 (a)- 
Parameters: r = 1, q 2 /[4(Ci + C 2 )] = 40I\ e d = 0, k B T = 5I\ 
g/j^BB = O.ir and capacitance asymmetry 77 = 0.5. Note that in 
our chiral system spin indices are included in the lead indices for the 
fluctuation relations. 



presence of arbitrary interactions. Promising avenues for fu- 
ture work include finite-frequency calculations and ac fields. 
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